Nothing moves a surface: vacancy mediated surface diffusion 
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We report scanning tunneling microscopy observations, which imply that all atoms in a close-packed 
copper surface move frequently, even at room temperature. Using a low density of embedded indium 
'tracer' atoms, we visualize the diffusive motion of surface atoms. Surprisingly, the indium atoms 
seem to make concerted, long jumps. Responsible for this motion is an ultra-low density of surface 
vacancies, diffusing rapidly within the surface. This interpretation is supported by a detailed analysis 
of the displacement distribution of the indium atoms, which reveals a shape characteristic for the 
vacancy mediated diffusion mechanism that we propose. 
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Mobility of close-packed metal surfaces is usually 
thought to be restricted to the direct vicinity of steps, 
where atoms most easily 'come and go'. The diffusion 
of adatoms along or between steps leads to characteristic 
step fluctuations. Also adatom and vacancy islands are 
known to move via rearrangements at their perimeter. 
Many STM-studies have been devoted to the mobility 
of surfaces. Most of these have focussed on the motion 
of steps [fy-[|, islands M or adsorbates M. Recently it 
has been suggested that also surface vacancies can be re- 
sponsible for self-diffusion of islands |6| . Unfortunately, 
there are no experimental techniques available with both 
the spatial and temporal resolution necessary to follow 
the diffusion of naturally occurring vacancies in a close- 
packed metal surface. 

Indium which is deposited on Cu(001) has been found 
to modify the epitaxial growth of copper on this sur- 
face. Rs presence results in layer-by-layer growth instead 
of rough three-dimensional growth J7J . After deposition 
the indium atoms proceed to steps on the copper surface 
[|S|,|9| . At temperatures just below room temperature they 
are incorporated in the outermost layer on substitutional 
terrace sites. In this study we have used indium atoms 
that are embedded within the first layer of a Cu(001) sur- 
face to monitor the diffusion of surface atoms |1(J . Our 
observations lead us to conclude that surface vacancies 
are responsible for the mobility of the indium and that 
this close-packed metal surface is far from static, even at 
room temperature. 

The experiments were performed with a variable tem- 
perature scanning tunneling microscope (STM) llT[] in 
ultra-high vacuum (UHV). A Cu single crystal of 99.999 
% purity was mechanically polished parallel to the (001)- 
plane |12|. Prior to mounting the crystal in the UHV 
system we heated it in an Ar/H2 atmosphere to remove 
sulfur impurities. The sample surface was further cleaned 
in UHV by several tens of cycles of sputtering with 600 eV 
Ar ions and annealing to 675 K. After approximately ev- 
ery fifth cycle the surface was exposed to a few Langmuir 



of O2 to remove carbon from the surface. STM images 
showed a well-ordered surface with terrace widths up to 
8000 A. Small quantities of clean indium were deposited 
on the surface from a Knudsen cell. 

The starting point of the observations is shown in 
Fig. nl At room temperature we have deposited 3% of a 
monolayer of indium on the Cu(001) surface. The STM 
image shows a region around an atomic step, separating 
two flat terraces of the copper surface. The image was 
taken 42 minutes after deposition and shows that most in- 
dium atoms are within 150 A from the step. From the ap- 
parent height of 0.4 A of the indium atoms, we infer that 
they are embedded within the first copper layer. What 
we know from lower-temperature STM experiments is 
that newly deposited indium atoms first 'hop' over the 
surface until they encounter a step. At the step they 'in- 
vade' the outermost copper layer (both on the upper and 
on the lower side of the step), after which they diffuse 
away from the step, whilst remaining embedded within 
the copper surface layer. 

We follow the diffusion of the embedded indium atoms 
in the copper terrace by making series of images of the 
same area on the copper surface to form an STM-movie 
of the motion [la ]. To our initial surprise, we found that 
the indium atoms move via long jumps of more than a 
single lattice spacing, separated by long time intervals 
]14| . In addition, the movies show that there is a strong 
tendency for nearby indium atoms to jump at the same 
time. Fig. illustrates this peculiar motion with a set 
of three images taken from a movie measured at 320 K. 
From the STM movies we have measured the distribution 
of jump lengths of the embedded indium atoms, which 
has been plotted in Fig. |[ Note that there is a significant 
probability for the indium atoms to make jumps as far 
as five lattice spacings. 

The long jumps and the high probability of nearby 
indium atoms to jump simultaneously, suggest strongly 
that diffusion of the indium is mediated by another par- 
ticle, which diffuses so rapidly that it remains invisible 




FIG. 1. A 548 x 409 A 2 STM-image of a step on a Cu(001) 
surface, taken 42 minutes after deposition of 0.03 monolayer 
of indium at room temperature. Embedded indium atoms 
show up as bright dots. The image shows a high density of 
embedded indium atoms near the step. (I t = 0.1 nA, V t = 
-0.70 V) 




FIG. 3. The distribution of jump vectors measured from 
STM-movies at 320 K. Plotted is the probability for jumps 
of an indium atom from its starting position to each of the 
shown non-equivalent lattice sites. To illustrate the unusual 
diffusion behavior, the expected jump vector distribution for 
the case of simple hopping is plotted to the right. 



to the STM. The scenario that we propose is that the in- 
dium moves over several lattice spacings during a multi- 
ple encounter with a single assisting particle by changing 
places several times with that particle. The two obvious 
candidates for this particle are adatoms (copper atoms 
on top of the surface layer) and vacancies (missing atoms 
in the outermost copper layer). We can rule out the first 
possibility on the basis of Fig. pi. If an indium atom were 
to change places with an adatom, it would thereby be- 
come an adatom itself. We know from Fig. u\ and from 
other observations that indium adatoms rapidly hop over 
the outermost copper layer to the steps, without entering 
the copper surface directly. This means that if an em- 
bedded indium atom would trade places with a copper 
adatom, it would immediately disappear from the STM 
image and reappear somewhere at the step, which is def- 
initely not what we observe. 

Fig. [| illustrates how a single surface vacancy can dis- 




FIG. 2. Three 50 x 50 A 2 STM-images selected from 
a movie measured at RT illustrating the unusual diffusion 
of embedded indium atoms. In the time interval of 160 s 
between images (a) and (b), no diffusion of the embedded 
indium atoms has taken place. In image (c), taken 20 s later, 
a diffusion event has taken place. Both indium atoms present 
in images (a) and (b) have moved over several lattice spacings 
and two more indium atoms have jumped into the imaged 
region, (I t = 0.9 nA, Vt = -0.58 V) 



place an atom in the outermost copper layer, either an 
indium or a copper atom, over several lattice spacings. 
In this mechanism, the length of the long jumps of the 
indium atoms depends on the average number of times 
that a single vacancy changes places with an indium 
atom, and we associate the frequency of the (long) in- 
dium jumps with the frequency with which the indium 
is encountered by new vacancies. We have measured the 
distribution of time-intervals between consecutive jumps 
(see Fig. ||). The waiting time distribution is purely ex- 
ponential, from which we infer that individual long jumps 
are uncorrelated in time and are therefore caused by dif- 
ferent vacancies, independently formed at random times. 
The fact that a single vacancy will usually encounter var- 
ious In atoms, naturally explains the tendency for nearby 
indium atoms to jump at the same time. 

The fact that we never see individual vacancies in the 
STM images and the fact that the STM movies do not re- 
solve the elementary steps in a multi-lattice-spacing jump 
need not surprise us. Using the Embedded Atom Model 
(EAM), we estimate that the formation energy of a va- 
cancy in the Cu(001) surface is 0.51 eV and that the ac- 




FIG. 4. A ball model illustrating the vacancy-mediated 
diffusion mechanism of the indium atoms. 
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FIG. 5. Time-interval statistics for subsequent jumps of 
individual indium atoms, measured from STM-movies at 320 
K with a time per image of 1.88 s. The dotted curve is an 
exponential fit with a time constant r = 8.5 s. 



tivation energy for a surface atom to exchange with the 
empty site, and thereby move the vacancy, amounts to 
0.29 eV |15[. Based on these estimates, we expect that at 
room temperature only one surface atom out of roughly 
6 • 10 9 is missing, and that each empty site changes posi- 
tion with a high frequency, on the order of 10 8 Hz. These 
numbers are typical for close-packed metal surfaces and 
illustrate why it is so difficult to sec the vacancy diffusion 
at all. At low temperatures, where vacancy motion would 
be slow enough to be followed by an inherently slow in- 
strument such as the STM, the probability of finding a 
vacancy is hopelessly close to zero. At temperatures high 
enough for the surface to contain a sufficiently high den- 
sity of vacancies, the vacancies move much too fast to be 
imaged at all. 

In order to obtain a quantitative understanding of the 
jump vector distribution of the embedded indium atoms, 
we performed a numerical calculation as well as a con- 
tinuum approximation, according to the following model: 
The Cu(001) surface is a finite I x I square lattice, with 
copper atoms at the lattice sites; the boundary of the 
lattice corresponds to the steps. One copper atom at the 
center of the lattice is replaced by indium, and a vacancy 
is released one atomic site next to it. The vacancy per- 
forms a biased random walk, its hopping probabilities to 
the four different directions from each site are set from 
the diffusion barriers calculated using the EAM poten- 
tials jn|. 

The vacancy displaces the atoms in its path, includ- 
ing the indium atom. When the vacancy arrives at the 
boundary of the lattice, it is annihilated (it recombincs 
at the steps). At this moment the displacement of the 
indium atom is evaluated, and the whole process is re- 
peated for the next vacancy to acquire the distribution 
of the In jump vectors. 

For the case of equal diffusion barriers and infinite lat- 
tices, this problem has been solved analytically m&. Al- 



though the results in some limits are quite similar to 
our continuum solution (see below), the equal-barrier re- 
sults are not directly applicable to the case of indium 
in copper. Instead of moving isotropically, the vacancy 
neighboring the indium atom has a much stronger pref- 
erence to jump towards the indium than to other direc- 
tions, based on EAM barriers at room temperature (see 
below). This difference has a significant impact on the 
indium jump distribution: the mean square displacement 
is about 2.2 times larger than in the equal-barrier case, 
while the overall shape of the distribution is about the 
same p~5| . 

It is computationally beneficial to separate the motion 
of the vacancy from that of the indium atom. For the 
indium atom, only the direction of the next return of the 
vacancy is of importance, rather than the vacancy path 
which leads to it. Therefore it is enough to calculate the 
probabilities of first return of the vacancy to the indium 
atom from the four different directions after it left the 
indium in one direction, as well as the probability of the 
vacancy's recombination before return. The In atom per- 
forms a random walk, where the direction of each step 
with respect to the previous one is chosen according to 
these return probabilities, and after each step the walk 
terminates with the vacancy's recombination probability. 
This procedure yields the proper final jump distribution, 
while giving up the time information, which is experimen- 
tally irrelevant anyway. (This approach is valid under the 
assumption that the environment of the indium does not 
change with the steps it takes, i.e., it is still close to the 
middle of the lattice.) 

In practice, instead of using Monte-Carlo type meth- 
ods, we enumerate the possible trajectories to obtain the 
return probabilities and the indium jump vectors; this 
provides superior convergence. The following numerical 
values were calculated for T = 320 K (EAM-barriers) 
and for a lattice size I = 401, which corresponds to 
the typical experimental terrace width of 1000 A. Af- 
ter leaving the indium atom to the right, the vacancy's 
return probabilities from the four directions are the fol- 
lowing: p right = 1 - 2.4 • 10~ 7 , p up = Pdown = 1.1 • 10~ 7 , 

Pleft =4.2- 10 , and the vacancy recombines with proba- 
bility p ICC = 1.1 • 10~ 8 . These values depend very weakly 
on the lattice size I. The fact that two dimensions is 
the marginal dimension for the return problem of a ran- 
dom walker implies a logarithmic / dependence of p tcc . 
The root mean square jump length of the In atoms is 3.5 
nearest neighbor spacings. The full distribution of the In 
jump lengths is plotted in Fig. ra together with the exper- 
imental values. The quantitative agreement supports our 
interpretation of the mechanism of the indium diffusion. 
Finally, we show that a simple continuum approach 
to this problem gives a quite good approximation to the 
jump statistics. Let us denote the probability of "mobile 
indium" at position r with g(r, n) , where the counter 
n measures the number of times the vacancy returns to 
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FIG. 6. The distribution of the jump lengths of the indium 
atoms at 320 K. Filled circles correspond to the experimental 
values, open circles are from the numerical calculation, and 
the solid curve is the isotropic continuum (Bessel) approxi- 
mation. For comparison, the dashed curve shows the best-fit 
Gaussian distribution. 



the In atom. The indium is considered "mobile" while 
the vacancy is still around, and "immobile" after the va- 
cancy has recombined. The effective diffusion equation 
for g(r, n) is 



de( Z> n) =DV>e-eQ. 

an 



(1) 



The first term corresponds to the vacancy mediated dif- 
fusion of the mobile indium (isotropic, in this continuum 
approach) , and the second term to the recombination of 
the vacancy, which makes the indium immobile. The 
solution in case of Dirac-delta initial conditions at the 
origin is 



e(r,ra) = 



1 



inDn 



(2) 



We are interested in the final, "immobile" distribution 
of In: 



p(r) 



eg (r,n)dn = ±-^K \ 



(3) 



where Ko is the modified Bessel function of order 0. The 
parameters can all be calculated: e is the recombination 
probability p vcc of the vacancy, and the effective diffusion 
coefficient D can be calculated from the return probabil- 
ities, as will be discussed in detail elsewhere Jig] . 

This analytical solution gives a good approximation of 
the jump length distribution, as shown on Fig. a We 
emphasize that neither the results of the numerical cal- 
culation, nor the continuum solution contains any fitting 
parameter: everything is calculated from the EAM bar- 
riers, the temperature and the average terrace width. 

In conclusion, the diffusive motion of the indium atoms 
can be explained by the presence of a low density of ex- 



tremely mobile vacancies in the first layer of the sur- 
face. This interpretation is supported by the shape of 
the distribution of measured jump lengths. The root 
mean square jump length can be reproduced accurately 
in calculations if we take into account the chemical dif- 
ference between the indium and copper atoms. The the- 
ory further shows that the multiple encounter of a single 
vacancy with a copper atom in a clean copper surface 
should result in a root mean square displacement of the 
atom of 1.6 nearest neighbor spacings. Combining this 
number with the observed average jump rate of the em- 
bedded indium atoms, we calculate a diffusion coefficient 
for copper atoms in a Cu(001) surface of 0.42 A 2 -s _1 . We 
see that close-packed terraces of metal surfaces, such as 
Cu(001), cannot be considered as static, even at room 
temperature. The naturally occurring vacancies lead to 
a continuous reshuffling of the surface, as if it were an 
atomic realization of a slide puzzle! 
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